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1. Introduction 

Applications of the AdS/CFT correspondence to non-static geometries continue to be an 
active area of research. This is motivated partly be the inherent interest in non-equilibrium 
processes in strongly coupled quantum theories and partly by the pressing questions arising 
in this context in connection with the phenomenology of heavy ion collisions ||l| (and more 
generally of nuclear matter in dynamical situations). 

An important connection between slowly evolving geometries on the gravity side and hydro- 
dynamic states of N=4 supersymmetric Yang-Mills theory was uncovered in |^ following 
the earlier observations of . The simplest case of this connection involves only the metric 
in the bulk, and describes the hydrodynamics of a fluid without any conserved currents 
beyond the energy-momentum tensor. This so-called fluid-gravity correspondence was sub- 
sequently generalized in various directions 0, ^, §|, • 

A very important case is that of hydrodynamics with conserved charges. Apart from purely 
theoretical interest, the motivation for this arises in the context of applications of relativistic 
hydrodynamics to the evolution of quark-gluon plasma created in heavy ion collisions. 
In that case the conserved charge is baryon number^. More generally, one is interested 
in the properties of gauge theory plasma at finite temperature and chemical potential. 
Such a system is dual to gravity interacting with a vector field. The generalization of 
fluid-gravity duality to this situation was taken up in 0, ^, where a solution of the dual 
gravitational theory was described (up to second order in the gradient expansion) and 
transport coefficients of the charge-bearing plasma were calculated. Since the bulk gravity 
theory as it appears in a string compactification]^ includes a Chern-Simons term, the dual 
hydrodynamics includes terms unexpected from the point of view of classical expositions 
of the subject. This has lead to very interesting developments which clarify the effects of 
anomalies in hydrodynamics!^. 

The major goal of this article was to determine the causal structure of the dual spacetime. 
This is motivated by two issues. First of all, the gravitational solution found in the gradient 
expansion is smooth apart from a singularity at r = (the black brane singularity). On 
general grounds it is expected that this singularity should shielded by an event horizon. 
Establishing this is important and nontrivial. The second issue is that hydrodynamic en- 

^In general one may also consider additional conserved charges, such as strangeness (as long as weak 
interactions are negligible). 
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tropy currents [10, 11, 12 1. The event horizon defines an entropy current and the condition 
that its divergence be non-negative imposes constraints on the transport coefficients^. It 
is natural to also consider entropy currents defined by dynamical horizons in the geometry 
1 15, 16, 1^]. This subject will be taken up elsewhere [p^ ]. 

Much of this article is devoted to establishing the dual solution itself. Compared to the 
earlier works the calculation presented here differs in a number of ways. First of all, 
the results generalize the findings of |Q, ^ in that an arbitrary weakly curved boundary 
metric is admitted. As in the case without charge, the result can not be obtained just by 
covariantizing the flat-boundary result, since there is an additional term which involves the 
curvature of the boundary metric Q . More importantly from the point of view of locating 
the event horizon, the form of the solution obtained here is somewhat simpler and more 
explicit than the earlier results. This is due to a number of factors. Since one of the goals 
of this study was to determine the causal structure of the dual spacetime geometry, it 
was convenient to use a parameterization of the zeroth-order solution such that the event 
horizon at that order can be expressed in a simple way. Such a parameterization (found in 
||T^) also simplifies the gradient expansion, and especially the perturbative determination 
of the event horizon of the slowly evolving geometry. Another technical difference is that 
the present paper makes use of a different gauge than that used in |^, |5|, following the 
choice made in [||, |^. The advantage of this gauge, apart from simplicity, is that ingoing 
null geodesies are simply curves of constant boundary coordinate x. This defines a natural 
bulk-boundary map, which is an important element of the holographic construction of the 
hydrodynamic entropy current. Finally, explicit Weyl covariance in maintained throughout 
(as in [^, 0]), which simplifies the calculations as well as the form of the final results. It 
turns out, that taken together these simplifications make it possible to write completely 
explicit formulae for the metric and gauge field, which are somewhat more complicated 
than those given in 0] for the case of uncharged plasma, but not significantly so. In 
particular, it is manifest that the results given here reduce to those of ||7| in the limit of 
vanishing charge. Following the standard procedure of holographic renormalization simple 
and complete expressions for the transport coefficients are also obtained. 

The structure of this article is as follows. Section |2| reviews the relevant static gravity 
solution and some aspects of its thermodynamics. Section ^ describes the general form of 
the solution and an overview of the computation. Section presents the solution up to 



2 See however ^ 



-3- 



second order in the gradient expansion. Section ^ describes the results of the holographic 
renormalization procedure and the results for the transport coefficients. The causal struc- 
ture of the geometry is studied in section |6|, where the event horizons are located. Some 
closing remarks follow in section |^. 

2. The bulk theory, black branes and thermodynamics 

The action of the five-dimensional Einstein-Maxwell theory under consideration reads 

S = ^J d'x^g [§ + R-\f' + ^6^^^^^^^FbcFde) (2.1) 

where the cosmological constant is denoted by For this theory to be a consistent 

truncation of type IIB supergravity^] the Chern-Simons coupling k has to assume the 
value I/2V3. 

The action leads to the equations of motion 

Gab - Q9AB + 2FacFE + ^qabFcdF^^ = 0, 

VbF^^ + Ke^'^^'^'^FncFDE = 0. (2.2) 

Static black hole solutions posessing spherical symmetry as well as their thermodynamics 
were discussed in detail in |^]. As shown there (following ||2^) a scaling limit gives rise to 
the following solution: 

At = h{r) , (2.3) 

where 

/«-(i-jl)0 + |-|i). M 

This solution possesses a planar event horizon. In consequence of AdS/CFT duality it 
describes the thermodynamics of plasma in equilibrium in flat Minkowski space. 

The position of the (outer) event horizon is at r = tq. There is also an inner horizon at 
r = r_, where 



1 



,2 



-- = ^-oh 1 + 4^-1 • (2-6) 
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Using standard Euclidean techniques one finds the Hawking temperature of the outer hori- 
zon m 

^l-^)- (2.7) 



7rL2 V 2ro 



Note that the temperature vanishes for the extremal black hole with q /r^ = 2. 



As discussed in [19| the chemical potential in the field theory is related to the asymptotic 



behaviour of the gauge field. It can be related to and the Hawking temperature by 



ro = vrL^- | 1 + ^/1 + ^^j . (2.8) 



For the purposes of fluid-gravity duality it is appropriate to use coordinates which are not 
singular at the event horizon. The Schwarzschild-like coordinates used in suffer from 
a coordinate singularity there. The choice made in was to use Eddington-Finkelstein 
coordinates. Starting from the coordinates used above one can use a transformation of the 
form r = r' + F{r') to reach such a gauge. In the present case this results in the following 
expression of the charged black brane solution: 

ds^ = 2drdt - ^-^^df + +^{dx^ + dy'^ + dz^) , 

A = ^, (2.9) 

2r^ 



where the function / is given above in eq. (2.4) 



3. General form of the solution 
3.1 The gradient expansion 

The method of ||2| mimics the way relativistic hydrodynamics arises from the static, ther- 
modynamic description. The energy-momentum tensor of perfect fluid hydrodynamics is 
just a boost of the equilibrium energy-momentum tensor, where the temperature and boost 
parameters are allowed to depend on position. Following the same idea one considers the 



boost of (2^), which describes equilibrium states, to some constant velocity n, and then 
allows this velocity and the temperature to depend on x. The boost parameter is a 4- 
component velocity vector in the directions, normalized so that u^u^^ = — 1 in the sense 
of the boundary metric h^y (metric on the conformal boundary of the locally asymptotically 
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AdS spacetime ( |3.lD ). Thus one is lead to consider the geometry 

ds2 = (p^^ _ 2Bu^Uu) dx^dx*" - 2?x^dx^dr, (3.1) 



where^ 

1 



and 



P^iu = + u^Uy (3.3) 

is the projector operator onto the space transverse to n^, The vector potential takes the 
form: 

A = ^n^dx^ (3.4) 

The constant parameter b is just l/rg in the notation of the previous section. The solution 
described there is recovered by going to the frame where {u^) = (1,0,0,0). 

The geometry described above has a curvature singularity at r = 0. The latter is shielded by 
the event horizon at r = 1/6. The parameter h appearing in ( p.ip is related to the Hawking 
temperature T of the event horizon by eq. ( p.?] ), which in the notation introduced above 
reads 

The lines of constant in ( |3.lD are ingoing null geodesic, for large r propagating in the 
direction set by u^, and the radial coordinate r parameterizes them in an affine way |ll| . 
Unlike black holes in asymptotically flat spacetime, the metric ( |3.1| ) supports perturbations 
varying much slower within the transverse planes than within the radial direction. The 
parameter controlling the scale of variations in the radial direction is b. 

The field configuration described above is a solution of the equations of motion for constant 
6, q, u^. If these parameters are allowed to depend on x, the equations are violated by terms 
proportional to gradients of 6, g, u. To cancel these, so as to ensure that the fields still satisfy 
Einstein equations corrections need to be added to the metric and gauge potential order 
by order in an expansion in the number of gradients. Thus, if 6, and h^y are allowed 
to vary slowly compared to the scale set by b, the metric ( |3.l| ) should be an approximate 
solution of nonlinear Einstein's equations with corrections organized in an expansion in the 
number of gradients in the x^ directions. As in the uncharged case |^], this turns out to be 

^From now on the constant L is set to unity, and the notation is chosen to resemble that of reference |7|. 
*The notation is chosen so that in uncharged hmit (g — > 0, k — >■ 0) B is equal to B{br) as defined in |7|. 
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possible if and only if b, q, u satisfy differential equations which can be interpreted as the 
equations of hydrodynamics |2|. 

Technically this can be done by considering an arbitrary point, say x = and expanding 
in a Taylor series 

n^'(x) = n^(0) + ex'^^cu^'{{)) + ... 
b{x) = b{0) + ex"da 6(0) + . . . 

q{x) = (7(0) + ex"5„ q{0) + ... (3.6) 

Each derivative with respect to the "boundary coordinates" x is tagged with a power of e 
for power counting purposes, and e is set to unity at the end of calculations. 

The four-dimensional boundary metric components h^i, are also expanded around x = 0, 
assuming that at zeroth order h^^{0) = rj^^. Moreover, to simplify computations it is very 
useful to adopt a locally geodesic coordinate system on the boundary, so that all first order 
derivatives of hf^i^ vanish. Thus one has 

In general, the second order derivatives cannot of course be set to zero in this way, but 
these contributions are guaranteed to be tensorial and so the final results obtained are 
covariant in the boundary sense. 



3.2 Weyl covariance 

Weyl covariance in the bulk arises as an extension of the conformal symmetry of N=4 
supersymmetric Yang-Mills theory and its charged cousins |10, |^]. A beautiful formalism 



allowing manifest Weyl covariance was introduced by Loganayagam |10| and applied to 
fluid-gravity duality in [^, |^] . 

The calculation is simplified considerably by adopting the approach of which imposes at 
the outset the conditions of Weyl invariance on the possible form of the solution. Conformal 
symmetry of the dual field theory can be extended to the bulk as follows ^: 

9^,u e-^^g^,y, ^ e<^M^, b e'H and r ^ e^r (3.8) 



where (j) depends on the coordinates x^ 10]. The leading order metric (3^) is Weyl-invariant 



but due to the presence of dr it does not retain its form at higher orders. It can however 
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be written in a manifestly Weyl-invariant form upon introducing a vector field Au defined 

by 

This quantity is of order one in the gradient expansion and transforms as a connection 
under Weyl-transformations 

Au^Au + duCp. (3.10) 



As explained in [10| this connection allows one to introduce a Weyl-covariant derivative T>. 



The gauge field A is a vector field of Weyl weight zero. It will be taken in the gauge Aj. = 0. 
The general Weyl-invariant form of the metric ( |3.lD reads 



ds^ = (r^P^^ - 2r^Bu^u^) dx^dx'' - 2u^dx^ (dr + rAudx") . (3.11) 

The correction involving A compensates an inhomogenous term in the transformation of 
dr. It is of first order in gradients; further contributions are needed for a complete solution 
at this order, but they are by themselves Weyl-invariant. 



The static metric (|3.1| ) is a leading order approximation to a spacetime whose metric is of 
the formQ 

ds^ = (g^^ - 2Uf,Vu) dx^'dx" - 2'u^dx^ (dr + rAudx") , (3.12) 

with the condition u^G/iu = completely fixing the gauge freedom This implies a choice 
of gauge which is different than that of fs], |^ . It is easy to see that lines of constant x are 
geodesies (affinely parameterized by r), as in the case of the static metric. 

The metric ( |3.12| ) is manifestly Weyl-invariant, provided that the functions are of unit 
Weyl weight and ^^i, are Weyl invariant. The simplest way to construct them is by summing 
individual tensorial contributions of appropriate Weyl weight order by order in the gradient 
expansion, multiplied by scalar functions of the Weyl-invariant combinations rb and qb"^. 

A powerful tool for generating Weyl-covariant gradient terms is the Weyl-covariant deriva- 
tive V^, which uses the connection Af^ (|3.9|) to compensate for derivatives of the Weyl 
factor coming from derivatives of Weyl-covariant tensors. It has the property that a Weyl- 
covariant derivative of a Weyl-covariant expression is itself Weyl-covariant with the same 
weight (see Appendix ^ or the original publications |10, |l|, |7| for details). 



-8- 



4. The solution order by order 



4.1 Zeroth order 



The solution at leading order is just the boosted charged black brane solution (3J). The 
metric is obviously of the form ( 3.12| ) with 



G^u =r'Pf,, (4.1) 

(strictly speaking, to retain only leading order terms, the term involving the Weyl connec- 
tion in ( p.l2 ) must also be dropped, as it is of first order). 



4.2 First order 

To find the complete solution at first order one needs to classify the possible terms that 
may appear. At first order there are: 

• no Weyl-covariant scalars 

• one Weyl-invariant pseudo vector: 

= Cf^^xpu'D^P (4.2) 

• one Weyl-invariant vector: 

= q-'PHT^uQ (4.3) 

• one Weyl-invariant symmetric tensor of weight w = —1: 

o-pu = T^(f,Uu) (4.4) 

Therefore the general form of the Weyl-invariant structures appearing in the metric is 

= r^Bu^ + rFil^ + br^FoVo^, 
G^,. = r^P^,. + 2br^F2a^,, (4.5) 

where Fi and F2 are functions of the Weyl-invariants rb and qb^, which are to be determined 
by solving the field equations up to linear order in gradients. The factors of b and r above 
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were chosen to ensure the correct Weyl weights (and partly also for convenience - the choice 
made above leads to simpler differential equations). 

Similarly, Weyl covariance implies that the vector potential takes the form: 



^ = [^^^^ + ^O^A* + ^O^OmJ da^''- (4-6) 

Once again, Yq is a function of rb and qb^ and will be determined by equations of motion. 

To find the unknown scalar functions ^1,^2,10 the Ansatz described above is inserted into 
the field equations and expand in gradients using (|3.6[) . This process, while tedious^, is not 
significantly more complicated than what has to be done to reproduce the results of 0. 

As discussed in the bulk field equations are of two types: the constraint equations, which 
impose consistency conditions on the slowly- varying background parameters b{x), q{x), u^{x) 
and the dynamical equations which determine the functions which appear in the field 



Ansatz (13.121) 



The constraint equations take the form 



d^q = -qdiUi, (2 - b'^q^)dib - 26(1 + b^q^)doUi = qb'^diq, 

dob = ^bdiU^. (4.7) 



one can rewrite this as 



In terms of the Weyl-covariant derivatives 

V^b = d^b - A^b, 

'Df.q = d^q + 3A^,q , (4.8) 

(b^'^ - 2)Pf"'V^b = qb'^P'^'V^q, 

u^'V^b = 0. (4.9) 

These equations can be interpreted as the equations of hydrodynamics at order zero, which 
describe the perfect fluid limit of the supersymmetric Yang- Mills plasma on the boundary. 
This will be discussed in more detail in section (^) below. 

The remaining differential equations can be solved for the functions Fq, Fi, F2, Yq, Yq- Tak- 
ing into consideration the requirements: regularity across horizons and normalizability of 



^Thus best relegated to symbolic manipulation software 
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the metric one can fix all constants of integration to find 
_^ + 2 (ft^g^ + 2) 

° ~ 2br^ Ab^r^ 4r6 (66^2 + 1) + 

H ^TTT-S / dx- 

Fi = 



266r6 (l+^)2(_56g2+^2 + ^4)2 



f>00 

Fo = / dx 



x(l + X + x^ 



6r 



;i + x)(-6V + a;2 + x4) ' 



3b^Kq^ 



2r2 (56g2 + 1) ' 
~ ^ ^/35g(2 + bV) 

° 8(1 + 66^2)^2 



V363g /-"odx /-"^ /(-6V(l + 2j/) + j/2(3 + 2j/ + j/2)) 

2 A, X^y, (l + y)2(_56g2 + y2 + y4)2 • ^ • "J 

The symbols b and (7 appearing above are understood as values at an arbitrary point x, 
not necessarily x = 0, which was just an irrelevant choice made to implement the gradient 
expansion. As stressed in the original papers on fluid-gravity duality the equations are 
ultralocal in x, which is the key feature which allows one to determine the solution. 



In the uncharged limit (g — )• 0, k — )• 0) one obtains 

^(9=0) 



3-1 



which is exactly the function F{br) appearing in the solution presented in 0. In this limit 
Yq and Fi vanish (as they should). 

4.3 Second order 

As in the previous section, the first step is to determine all the relevant Weyl-covariant 
structures which may appear at second order. In doing this one has to discard any terms 
which vanish or are not linearly independent after the first order constraints ( |4.9| ) are taken 
into account. 

The results of this analysis are as follows: 
• Scalars: 

Si = b^a.^a'^'^, 
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52 = h'LO^.U^^ 

53 = b^n, 

54 = b^-^P^^^Di.qD.q, 

55 = b''q-'P^"'D^D,q, 
Se = b^q-^P'^^'l^D^q 

• Vectors: 





= bP^^Dpa^P, 




= hP^^DpU^'P, 


y^. 






= bq-^a^^'Daq, 


y^t. 


= bq~'^u^Daq 



(4.12) 

• Tensors: 

T21J.U = C^aul3lJ'°'lJ'^, 
Tzp,u = ^il^'^Xy +'^v^'y\fj., 
T4fiu = (^i^f^Xv - -^PixvC^afif^'^'' , 

Tbnu = ^fl^lJ^Xv + -^Pixv^a^lJj'^^ , 
^liiv — 2 \[p^ afiu)a 'J' " > 

Tg^, = q-^W^^^D^qD^q, 
Tg^, = q-^W^^D^D^q, 

Tiopu = q-^Ufj^Df^q, 

Here Il^f! is the projector which can be used to create symmetric, traceless tensors: 

= \ [p^pS + P^pS - Ip'^^P^..^) , (4.13) 
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and the parenthesis () denote symmetrization: 



A(^^^y.= A^^ + A^^. (4.14) 

On the basis of the above results one can write down the most general form for the fields 
allowed by Weyl invar iance: 



i=l i=l 
6 11 



The vector potential takes the form: 

^ = ( + ^o^A* + ^o^Om + ^ E ^i-^i"/^ + E dx'^ (4-16) 

V ^ i=l i=l J 

The 39 coefficient functions Ki,Li,Ni {i = 1,...,6), Wi,Yi {i = 1,...,5), and Hi {i = 
1, 11) all depend on the Weyl invariant variables b^q and br. Solving the equations of 
motion up to the second order determines these functions uniquely once regularity condi- 
tions and a frame choice are imposed. Moreover, from the constraint equations one obtains 
relations involving second-order terms in u^, b and q. As in other similar contexts, these 
coincide with the equations of hydrodynamics at order one. 

This procedure is more cumbersome than at first order, but as before one can simplify it by 
first establishing the constraints and the using them to simplify the remaining equations. 

In this way all the functions appearing in ( |3.12| ) are determined. The integration constants, 
as at first order, are all fixed by conditions of regularity, apart from the constants appearing 
in functions Ki, ...,Kq, Yi, ...,Y^, which can be fixed by choice of frame. The results are 
listed in appendix |A[ 



5. Holography 

The holographic dictionary of the AdS/CFT correspondence provides a a straightforward 
prescription for calculating physical quantities in the boundary theory. In practice the 
computation involves subtraction of divergences, which is done in a systematic way using 
holographic renormalization , 22]. 
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The expectation value of the energy momentum tensor is given by 

T^v = -^4^ hm {K^^ - KH^^ + 3F^^ - E^^) , (5.1 



8ttG5 



r—^oo 



where Hf^^y is induced metric at the the surface r = const, Hf^u is the extrinsic curvature and 
- the corresponding Einstein tensor. Carrying out the calculation requires a careful 
evaluation of the asymptotics of the integrals given in appendix One then finds the 
following result 

^""le^GsV 64 + ^3 + 52 + 52 + 52 + 

, 2r4^. ^ 46V(-1 + 6V(12^^-1)) ^ 2^{36V^ cg^ ^ 



^2^9^!/ -t- ^2 



+ i^Tg^v + -lTTlQ^,y ) , (5.2) 



where 

5V + 1 r dx 



ci 



2 ii -b^q^ + X + x2 ' 

/■oo 

cs = I dxps{x), cg = / dxp9(x), cio = / dxpio(x), 



00 



with ^8(3^)5 ^9(2^); Pio{^) are defined in the appendix |A| (together with the second order 
solutions) . 

The result ( ^.21 ) agrees with (for d = 4) in the limit g — t- 0, k — )• 0. From it one can read 
off the values of the transport coefficients. 



This formula ( |5.2D shows that the equations of state are p = 3e (as expected for a conformal 
theory) and 

^ '1 T,2„2 



U + bV) (5.3) 



IGttGn b^ 

(which is also consistent with conformal symmetry). This formula shows an interesting 
duality property, that is, invariance under the substitution 

1 

e ^ ^. (5.4) 
Under this transformation the Weyl invariant quantity qh^ is inverted: 

^ ^- (5-5) 
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This duality can be easily expressed in terms of T and ^. 

Holographic renormalization of the current proceeds according to the formula 

= hm (5.6) 

Inserting the solutions one finds: 

" SvrGsV 2 ^2(1 + 6V) 8(1 + bV)^^ °^ 8(1 + 6V)^ 



M^V^T/ 364^g2 2a4 (b V + l) + ^/36 V ^ 

(1 + 6VP 2(6V + 1)' ''^^ 1662 (&V + 1)' 



as (fe6g2 + 1) + (24^2 - 1) g3 - V3b3g ^^ 

The constants 04 and 05 are also defined in (0): 

2V^b^q {b^^q^ (3 - 2x3) + 56,2^2 (.3^2 + 2x - S) + 4^3) 



00 



a4 := / dx 



x2 (66,2 _ 2) (56,2 _ ^2 (3,2 + 1)) 

VSb^q (b^^q"^ - ^b^q^ + 2) SV^b^qFo ib^q'^ (Sx^ - 9) + x^ (Sx^ + 5) ) 
"5 := n , ^ + / dx ^ ^ 2 — ^ ^ 



266g2 _^ 2 

This formula contains the values of charge transport coefficients. 



6. Causal structure of the dual geometry 



This section is devoted to locating the event horizon in the geometry discussed above. It 
is important to verify that the metric singularity at r = is shielded by an event horizon, 
as expected on general grounds. As explained in |jll[ (see also |jl^, |l^) is possible to find 
the event horizon in the gradient expansion by taking advantage of the fact that at zeroth 
order its location is known. 

To determine the event horizon it is simplest to assume that it can be presented as the 
level set of a scalar function S{r^x). This function must be Weyl-invariant and can be 
written in the gradient expansion in terms of all the independent scalars available up to 
some order. Since the event horizon at order zero is at r = 1/6, one has 

S'(r, x) = 6(x)r — g(x) (6-1) 

where g{x) is a Weyl-invariant function expanded in gradients of 6, q: 

9{x) = go{x) + gi{x) + g2{x) + . . . (6.2) 
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Here denotes a linear combination of all Weyl-invariant scalars at order k in the gradient 
expansion. Thus go{x) is a function of Weyl weight zero, i.e. it depends only on the 
combination qb^. There are no Weyl-invariant scalars a order 1, and 6 at order 2, so one 
expects to find 

go{x) = X{qb') 
gi{x) = 

6 

g2{x) = Y,hk{qb^)Sk, (6.3) 

k=l 

where the Si are the 6 independent Weyl-invariant scalars ( 4.12| ). The functions hi and A 



will be determined in due course. Once this is done, the expression for the position of the 
event horizon will take the form 



I + • (6.4) 



The normal covector to a surface of the form ( |6.1[ ) is 

m = dS , (6.5) 

which up to second order in the gradient expansion is 

m = r db + b dr + dX. (6.6) 

It is convenient to write the normal in terms of the Weyl-covariant derivatives (given in 
(|4|), (U)). One then has^ 

m = {rV^b + X' {b^V^q + Sqb'^V^.b) )dx^' + b {dr + rAf.dx^') . (6.7) 

Since the event horizon is a null surface this normal must satisfy m? = 0. At leading order 
this condition has two solutions: A = 1, which is the outer event horizon at order zero, and 
A = r_, where r_ is given in ( |2.6D , which is the inner horizon of the charged black brane. 
At order two the condition fixes the constants hk appearing in (|6.4|). Carrying out this 
calculation with the metric obtained earlier one finds the outer event horizon at 

_1 3 {b^q^ - 2) Ki {b^q, l) + 1 (ft^g^ + 1)^X2 {b^q, l) + 3b^^K^q^ 

" 6 + 3b (b^q^ - 2f ' ^ b{b^^q^ - 3b^q^ - 2) ^ 

1 c fKi{b^q,l) b^^q^ {m^^q^ + 28b^q^ + 20)\ 



126(2 -6V) V6(6V-2) 32 (6V - 2)^ (66g2 + 1)2 



^Here A' = X'{qb^). 
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^ K,{b\l) bV(6V + 2) N 

U (66g2 -2) 4 _ 3512^4 + 4) j + 

4 (6iSg6 _ 3^6 ^2 _ 2) (53^^ 1^ ^ V36Vg3 (56^2 ^ 4) (^356^2 ^2) 

+ ;5 Sr. (6.8) 

46 (-512^4 + ^,6^2 + 2)2 ^ ^ 

The inner horizon at second order is given by 

r = r_ + hiSi + /i2S'2 + h^Ss + h^Si + /isSs + hf^SQ, (6.9) 

where the coefficient functions hi can be found in appendix |C|. As seen from the expressions 
given there, these functions (given in the form of integrals) diverge at r = r_. This is not 
surprising, since the corrections to the metric and gauge field are divergent there already 
at first order, indicating a breakdown of the gradient expansion. 

It may also seem interesting to consider the extremal limit, in which r_ = 1/6 so the inner 
and outer event horizons at zeroth order coincide. This corresponds to qb^ = \f2. The 
gradient corrections to this horizon include divergent integrals, so one must conclude that 
the gradient construction breaks down in this extremal limit [S]. 



7. Conclusions 



The main results of this article are the determination of the geometry dual to hydro- 
dynamics with a conserved current in an arbitrary weakly curved background and the 
determination of its causal structure. The calculations of the bulk field theory solution in 
the gradient expansion described here follow the pattern of earlier calculations of this type 
and have benefited from a number of technical insights accumulated in this field. This 
made it possible to present the dual description in a form suitable for studying its causal 
structure following the method described in [^] (see also fl^). The location of the event 
horizon is found explicitly in section (^, thus showing that the locus r = is not a naked 
singularity. As is well known, the event horizon can also be used to define a hydrodynamic 
entropy current for the dual field theory [[III]. 

Apart from the event horizon, the charged black brane geometry possesses also an inner 
horizon. The gradient corrections to the metric (as well as the gauge field) are singular 
at the location of the zeroth order inner horizon, indicating a breakdown of the gradient 
expansion. It is straightforward to formally determine the location of the inner horizon 
in the gradient expansion. The double integrals appearing in the formal expression for its 
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location are divergent. This is reminiscent of the remarks in ||5| concerning the extremal 
limit, in which the inner and outer horizons coincide. As noted in [|5|, the solution obtained 
in the gradient expansion develops singularities as the temperature T vanishes with q 
nonzero. This is apparently an indication that this limit is not described by hydrodynamics 
on the field theory side of the duality 24 1. It could however be that the extremal limit 



could be taken in such a way that an alternative gradient expansion would be valid. 

An obvious extension of this work would be to determine the entropy current associated 
with the event horizon. This is left to a future work [^], where the entropy currents 
associated with dynamical horizons will also be discussed. 

Finally, it would also be interesting to include background fields in this calculation, along 
the lines of |2^, where some important partial results are obtained. 
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A. Second order functions 

A.l Scalar sector 

h'^x V^x^ \ (1 + X + x^f 



Ni = q dx 



oo 



br 



2 



2V3r 2r3 ^/^.hr"^ ^1 + xf{-Wq^ + x'^ + x^f 
SVSb^K'^q^ V3q V^b'^K^q^ 



5r-9 (56g2 1)2 452^5 ^7 (56^2 _^ i) ' 

N3 = 0, 

_ [°° ( VSb^Foq^x"^ {b^q^(2x - 3) + 4x - 3) (x - brf{br + 2x) 

^ ~ kr r\x - l)2(x + 1)2 (66g2 + (_^2^6 + 3,4 + ^2 f + 

2^/3^0^9x3 (66g2 _ 2) + x2) (x - br f{br + 2x) 
^ 6r4(x- l)2(x + l)2(-g256 + 2;4 + x2)2 ^ 



b^q^x{x — brY{br + 2x) 



8\/3r4(x - l)2(x + 1)3 (66^2 + 1)2 ^_^2^6 + 3,4 + 3-2)3 



( 
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+ 2{x - ifx^ (2x^ + 18x^ + 21x2 + 18a; + 7) + h^^q^ {-Ax^ + 2x^ + 9x^ + x- 14) + 
+ h^'^q^ (4x^ + lOx^ - ITx^ - - 15x^ + 15x^ + SSx^ + llx - 28) + 
+ b^q^ (8x^ + 32a; - 49x^ - 21a;^ - 33a;^ + Ibx^ + Ux^ + 16a; - 14) )) , 
ga; {h^q^{2x + 1) - a;^ (x^ + 2x + 3)) (x - 6r)2(6r + 2x) 



f oo 

A^5 = / dx 



4V36r4(x + l)2(-g2b6 + ;;c4 + x2)2 

/ 245^^,K,H.r-/>r)^(5r + 2.r) 

+ ^ (524 8 (^5 + 2^4 ^ 3^2 ^ 24x + 18) + 

2r\x-l)x\x + lf{b^q^ + lf{-q^h^ + x'^ + x'^fy ^ ^ 
+ ft^g^x^ (-gx'^ - 2x^ - 19x^ - 9x^ + llx^ + 7x^ + 9x^ + ix^ + 4x + 5) 

- 2x^ (3x^ + 2x^ + x^ - 3x2 _ 2x - l) 

+ 6^^g^ (2x9 - 4x^ - 7x^ - 18x^ - 48x^ - 45x'' - 54x^ - 48x2 + 12x + 18) + 
+ 6^2g4^2 _ 3^11 _^ 2x^0 - Sx^ + 15x^ + 43x^ + 47x^ + 60x^ + 57x^ + 15x^ + 
+ 6x2 _ _ 27 



_ 1 _ 466^2^^ 
2 " 36V " 5r6 (66q2 + i)2' 

^3 = 0, 



/•oc 



^~ ' V " 3(X - 1)(X + 1) (66 g2 _ 2) (66^2 _ ^2 (^2 + l)) + 



4Fo2x5 (6V (2x2 - 3) + 2x2) ^^^4 

^ 3(^2 _ 1)2 (_^,6g2 + ^4 + ^2)2 + 3(^2 _ l)2(_56g2 + ^2 + ^4) 1, ^ ^ + 



+ 



6^^g^(7x2 + 6x - 15) + 6^2^4(^.6 ^ ^2x2 - 5) + 6^92(33.6 ^ 3^2 _ 243. + jq)) + 



+ b^fx^ (23,3 _ 3^2 _ 

12(x2-l)2(66g2_2)(56g2 + l)2(^6^2_^2(^2 + i))2V V 7 

- 26^2^4 _^ 3^2 _ 22^: + 16) + 26^g2 (^^7 _ 3^3 _^ g^2 _^ 2x - 8) + 4x (x^ - x2 + 6x - 6) ) 
2x^iV4" lOx^Ni' 2x'^Ni^ 



3x/363g sVSb^ x/363g/' 
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/•°° / 2x^Fo (&V(a;-2) + 2(x- 1)) 

iftr V 3 (x2 - 1) (-^2^6 + a.4 + ^2) + g (^2 _ 1) (56^2 + (56^2 _ ^2 ^^2 + 1)) + 



-I 2_ _| 2_ _| 2. 

SV^Pq SV^Pq VSPq 



00 



' " L ^V(x - l)x(a; + 1) (66g2 + i) (56^2 _ ^2 (^2 + " ^^/^ " 3^53^ " ^f^^ + 

6^ Kg 



^/3(X - l)x2(x + 1)2 (66g2 + 1)2 ^^6^2 „ ^2 (^2 + i))2 

+ (x^ + x^ + x"^ - x2 - X - 1) + h^^q^ (2x^ + 2x3 + 4x2 - 7x - 13) + 

- 6^2g4 ^2a:8 + 6x^ + 12x^ - Tx^ - ISx'' - ISx^ - 21x2 + x + 13) + 

+ iPq^x^ (4x^ + 4x^ + 4x^ - 2x^ - lOx^ - 22x'^ - Sx^ + 7x2 + + j j 

J. ^ ^ , J_ ( - , a;(l+x + 2x2-bV(l + x)) 

^ 2b^r^^b^r^J^ \ 6(l + x)(-6V + a;2 + x4) ^ 

2^''?' /Z./A2 x2-3x6 + 6V(a;'-6) 



x^ 



1 666^2^6 5^2 752^2^4 

^2 - - 5^12 (^6g2 + - Y^6V8 ~ 5rio (66g2 + 1) + 

^ bb^^q^ + 36^2 (4;«2 + 5) ^4 + 15^6^2 + 5 



if:. = 



2066r6 (66g2 + 1)^ 
1 



1262r2 ' 



i^4 = -Tr-TTTTT^r— TT + / ^X 



3r(46V + 4) V48 (^2 - 1)2 (56^2 + i)2 ^^6^2 _ 3.2 (3.2 + i))2 

+ b^^q^ (8x^ + 3x^ - 20x3 - 2x'^ + 12) + 

- 4x2 (2x^ + 15x^ - 18x^ +x^- 2x3 - Sx^ + 6x - l) + 

+ 6^2^4 ^_g^9 _ 15^8 ^ 3g^7 _ 4^6 ^ 24x^ + 15x^ - 64x3 - 12x2 + 24) + 

- 46^g2 (4^9 + 153.8 _ 273.7 ^ 2x^ - 6x^ - 6x^ + 17x3 + Sx^ - 3) ) + 

+ 5 ^ ^ ( b'\^ (3x2 + Sx-7) + 

6 (x2 - 1)2 (66^2 + 1) (_^256 + 3.4 + 2.2)2 V V ; 

+ 6^2g4 (_i2;;c6 _ + 21x^ + 9x2 + 93, _ 14) + 3 (3^10 _ 4^6 ^ ^2^ ^ 
+ 6^g2 (93,10 _ 24a.6 _ ig^s ^ 21a;'' + 9x2 + 6x - 7) ) + 
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Fo^x^ (b^'^q^ (2x2 - 3) + h^q^ {-I2x^ + Ihx'^ + - 3) + 2 (2,x^^ - Qx^ + x^)) 
^ 3(x2 - if {-q%^ + + x'^f ^ 

x^ {b^q^ - 3x^ + 1) N^' X {b^q^ (Sx^ + 6) - ISx^ + Sx^) iV4' 
3^/363g 3v^Pg 
(66g2 (a;2 ^ 2) - 3x6 + ^2^ 

^/363g 

_ b^q^ + 2 1 /"^ /x(x2(3x3 + 3x2 + 3x + l) -fe6q2(^^2^2x + 2)) 

" 3M4fcVT4) ^ 7b, V 6(x + 1) (-g266 + ^4 + ^2) + 

x^ (66g2 _ 3^4 ^ ATg" X (66g2 (52;2 + e) - ISx^ + 5x2) 

3^/363g 3V363g 
iV5(6V (x2 + 2) -3x6 + x2) 



V363g 



262^9 1 /"^ / -b^Kq 

Kg = /TT n . -■^ + HIT / ^X 



V3r(6V + l) ' b^r^Jbr ^'^ ^2VS{x - l)x^{x + \)^Wq^ + if {b<^q^ - x^ {x^ + l)f 
■ {h^'^q^ (2x^ + 2x3 + x2 - 4x - 7) + 4x^ (Sx'^ + x^ + x^ + x^ - Sx^ - x2 - x - l) + 
+ b^^q^ (-8x^ - 6x6 ^ ;^Q^5 _^ 35^4 _^ ^^^3 ^ ;^q^2 _ 53. _ + 
+ 66^2^2 ^jg^io ^ 2x'^ + 32x^ + 5x^ - SSx^ - llx^ - 23x'' - Tx^ + 7x2 + 7^, + 5^ + 
+ 6^2g4j^Q^i2 _ 23.11 + 13^10 _ ;^4^9 _ 4g^8 _ 22x'^ + 

- 40x6 - 8x^ + 40x^ + 19x3 + 14x2 - x - 7) ) 

2V3b^Kq^Fo {b^q"^ - 5x^ + l) Nq (b^q^ {x^ + 2) - 3x6 + ^.2^ 

^ X (x2 - 1) (66g2 + 1) (66^2 _ ^2 (3,2 + 1)) + 



x^ (66g2 _ 3^4 ^ iVg" X (66g2 (52;2 + 6) - 15x6 + 5x2) jYg/ 
3V363g 3V363g 

A. 2 Vector sector 



) 



Jbr v/ a; 



_ /(2 + j/) 3^/36g 

'A. x3y, ''^i+y)2( ■ ■ + 

3^/36^-2.3 



(l + y)2(-56g2 + y2+y4)2 8(66g2 + i)r2' 



(66g2 + l)2^2' 

366^^2 r^dx /"^ ^3 (56^2 ^ 2) + 3y'' + 6?/^ + 9y^ + 6^/2 + 4^ + 2) 

¥^TTjbr ^ (y+l)2(-g256 + y4+y2)2 



Sb'^Kq'^ 



2r2 (66g2 + 1)2 ' 
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f°°dx f'^ 

Jbr J X 



y[ F , 2V3b^q{z - 1) 

' (1 - 2/2)2 (_56 ^2 + y2 + y4) ^^2 (^6^2 _ 3) (56^2 _ ^2 (^2 + 1)) 



^4 = - / -T dy- „.o,o, ,.0 , „,o , / dz 



+ 



Y5 = 



36^2g^ (z'^ + z + l) - h^q^z^ {z^ + z^ + 2z'^ -z + 3) + 2z^ + + 2z + 2) ) + 
20,4 + 1) + V3b\^ 

^36^9 /"°^dx /-""^ fydz/ 



6i2g4 (_2^2 + _ 10) + 2^2 (^4 - 1) + 26^ g2 ^^6 _ + _ 5) + 

+ 3zFo [h'^q^ (5^2 - 9) + b^^ (Sz^ + 10^2 - 9) + ^2 (Sz^ + 5)) ) + 

V3a5 + l) + St" (24k2 - l) _ 353^ 
^ 8V362r2 (66g2 + 1)2 

where 



2^36^^ (6^2^4 ^3 _ 22;3) + 56^2^2 /_3^2 _^ 2x - 3) + 4x^) 
04 := / dx 



poo 



X2 (66g2 _ 2) (66^2 _ 3.2 (3.2 + 1)) 



V363g (6i2g4 _ 356^2 ^2) /-oo 3 V363gFo {h^q^ (5^2 - 9) + a;2 (3x^ + 5) ) 
26V + 2 



7br 



(X + 1) (-^256 + 3,2 + ^4) 



2 (66g2 (2x2 - 3) + 2x2) r'^ y6 (2 + y) 



7, ^(l+y)2(_66^2 + y2 + y4)2; + 



3 (66g2^2 _ ^4^2 + ^2) 



W2 = -:r^ + 



453^5 (56^2 + 1) ' 

^2^4 



1 66V2''4 



26r r5(6V + l)^' 

.00 (;,6g2^2 _ ^4 /^2 ^ + ^2 Wg// (6Vr2 - 36"^ (g2 - r6) + r2) yg' 

W3 = 6r / dx[ - ^— — h -^^ 7^- — h 

Jbr ^ V^b'^qr^ V^b^r* 

VSnq {b (6r2 {Sb'^q'^ + 2) + r {Sb'^q'^ + 2) + Ib^q"^ - ?,b^r^ - 36^^ - 36^^ - 62^3) + 2) 



62r3 (66g2 + 1) {pr + 1) (64^2 _ ^,2^4 _ ^2) 



+ 



^2j,5 y ' 
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/■oo (56^2^2 _ 1^2 _ ^6) + ^2) {pr - 1) (6r + 1) ib\^ - b'^r^ - r^) Y^' 

W4 = br dx[ 7^- — — h 

Jbr ^ VSb^qr"^ V3b*qr^ 

2b'^r^FQ 2r3 (68^2 + 52^ _ ^ (^^6^2 ^ 2) 



^,6g2^2 _ 54 (^2 + ^6) + ^2 26 (66q2 + l) (5^ _ J) (5^ + J) (^4^2 _ 52^4 _ ^2 

[,2^4 _ vJI (U^„2J1 _ niA /„2 J,\ , ^2\ v ' 



roc (6V-l)(6V-fe 
W5 = br dx( - ^ ^ . 

+ 263.5 (56,2 + 1) (,2,2^1) (54,2 2,4 2) {^^^W [h^^ " + T^) + r^) + 

+ (-2 (t^gV + r)' - 365gV {b\^ + 2) + 26V6 (fe^g^ ^ -^^ ^ g^4^2 (^6^2 + l)) ) + 



58,2,4 _ 56,2 (g2 + ,6) + 52 

A. 3 Tensor sector 



,4 ; 



/^i = -(brf 



/"OO X 1 /*x \ 



/•OO 

H2 = -2{brf / dx- 

Jbr ~ 



X 

-6V + a;2 + a;4' 



/•OO 



+ 12b'' kY {h\' (2 + x^ + x') + 3) + - 1) {2x' + l) 



2x/36llKgV2 roo x'^^x+l 
Ha = — : TTT-^ — / ax- 



(1- 




-6*5^2 + a;2 

x 




(1- 


X2)(- 


-b^q' + x2 

a; 




(1- 


X2)(- 


-66g2 + a;2 


+ x^) 



1 + x2(x + l)(-6V+a;2 + x4)' 

^^7 = 0, 

|-oo ^ /.X 

= -(^'r-)^ / dXy:; „2A/ I.fi^2 , ^2 , „4N / '^V P&iv) ^ 

Ibr 

H, = -{brfj^ d^ (l_,2)(_56,2 + ,2 + ,4) / 

dypio(y), 

i^n = 0, (A.l) 
where: 
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_ _4qx'^ {h^q^ + 1) {h^q^ (x^ - 3) + Zx^ + x'^) dFo 

^^^"^^ ''~ (x2 - 1) (66g2 _ 2) (66g2 _ (^2 + l)) + 

(x2-l)2(-g256+^4 + ^2)2 + 



,.."8 1 o^4 

(x2 - 1)2 (66g2 _ 2) (66^2 + 1) (_^256 + 3.4 + ^2)2 



+ ^^^g^ (2x'' + 9x^ - 22x2 - 12x + 24) + Ab^q'^x'^ (4x^ + 3x^ - Gx'' + 2x'^ - 9x + 4) 
- 6x^2 + b^^q^ (8x® - 3x^ - 6x^ + 8x^ + 9x^ - 28x2 + 6) + 
+ 26^2^'' (3x^2 + lox^ - 15x^ + 6x^ - 18x^ + Sx^ + 24x - 9) ) 



b^q^x 



4 (x2 - 1)2 (66g2 _ 2) (56^2 + 1)2 (56^2 _ ^2 (,^2 + l))2 

62^g8 (3x'^ - 12x3 + 8x2 _^ _ _^ g^2 ^3^6 _ q^5 + 3.4 _ 3^.2 ^ g^. _ 1) + 

+ b^^q^ (-3x^ - 36x'^ + 44x^ + 9x^ - 72x3 + 28x2 + 132x - 96) + 
+ 46^g2 (33.8 _ 24^7 ^ 15^6 _ 9^4 ^ 12x3 + 23x2 - 12x - 6) + 

- 66^2^4 ^^8 ^ 14^7 _ ig^6 + ^ iQ^3 _ 14^2 _ + ^4) ^ ^ 

^= (.2_i)(,6,2 + l)(,a,2_,2(,2 + l)) (^^^^ (^^^^(3- - 4) + 6x - 4) + 



+ 2xFo (6^2g4 (^^2 _ 3) + 56^2 (^3^6 ^ 2x2 - 3) + 3x^ + 
-6^/36^«;g3 



pio(x) := „ ^ x2 (2x5 - 3x4 _ 2;;c + 3) + 

x4(x2-l)(6V + l) (-«266+x4 + x2)V ^ ^ 
+ 6^g2 (^7 _ 3^6 _ 3^3 ^ g^2 + _ 7) _ ^12^4 ^^3 _ 3^2 _ 4^ ^ 7) ^ 

+ 4xFo (512^4 (^3^2 _ 4) _ ^6^2 (^^6 _ g^2 + 4) _ ^2 (^^4 _ 3) 



The existence of terms 1/(1 — x2) m. H\, H3, H4 may suggest that there is a problem 
at the outer horizon. However, examining the near-horizon behaviour one can find, that 
each of these functions is regular. To see this however, one must make use of the explicit 
form of F2. Just as at first order, all constants of integration were fixed by regularity, 
normalizability and choice of Landau frame. 
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B. Weyl weights 

Weight 4 : tensors from T|"^ to T^i , 

Weight 3 : q,a^"', w^'^ 

Weight 2 : g^'^J^', Vq'', vectors from Vf- to V^^, 

Weight 1 : T, /u, u^, r, 

Weight : Z^, Vq^, Weyl-invariant scalars, all covariant vectors and tensors 

Weight -1 : b,u^, a^y, uo^v 

Weight -2 : g^y. 



C. Inner horizon at second order 

The location of the inner horizon is given by ( |6.9D . The coefficient functions appearing 
there are given by 

/ii = 



254ri + 62r2 + 1 3 (52^2 ^ 2) (-26^1 + 62^2 + 1) 2 ' 

3 



y^rlK2 (^6r_ Y'&2r2 + 1, fer_ j 3^2 ^^2^2 _^ -|^^ 



-264ri + 62^2 + 1 ^_ (264^4 _ 52^2 _ (^55^4 _^ 53^2 _^ ^ 

^2^3 



12 (-264ri + 62^2 ^ 



7-3 ( ^3^-2 _|_ 5") 2 , 

''4-4 (52^2 + 2) 2 (64^4 + 52^2 + 1) (264ri - 62^2 _^y\ ^- + ^- ^« ^- + 

+ 3262r2 + 8 + 86^3^13 _^ 32512^12 ^ ig^n^n + 806^°ri° + }?r^_ + 1306V^ - llft^rl + 

54^5 ^^2 Uj.^Jh'^r'i + 1, 6r_) (br^Jpr'i + 1, 6r_ 

' -464ri + 262r2 +2 -264ri + 52^2 + ^ 



b^ri {Ab^r^_ + 76^ri + 562r2 + 2) Fq (hr^^h\'i + 1, 



6r_ 



(62r2 + 2) (-264ri + 52^2 ^ 



2 



^ (^6r_ ^62^2 + 1^ j ^3^4 (4510^10 + 11^8^8 + 2066r6 + 2164ri + 1262r2 + 4) 

^5 = -264rl + 62r2 + 1 + 4 (-264ri + 52^2 ^ 2 (56^6 ^ 3^4^4 ^ 352^2 ^ 2) 
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" -26Vi + 6Vi + 1 ^ 

_ V3Acr_ (62r2 + l) 3/2 (264ri + fcVi + l) [Ah^rt + 36^?. + 2) 
(-264rl + 62r2 + l) 2 (^6^6 + 3^,4^4 ^ 3^,2^2 + 2) 

V36Kr2 (62r2 + l) 3/2^0 (br^^jprl + l,6r_) 
268^8 + 66^6 - 262r2 _ 1 " 

In the last formula above, some functions, such as for example K4 (^br-^b'^r^ + 1, brJ^ are 
singular. In fact at the inner horizon, starting from first order, metric and vector potential 
are also singular. For instance, the first order function F2 can be rewritten as follows: 



(0.1) 



/ 



F2= dx 



br (l + x)(-56g2 + a;2+x4) 

hr (l + x)(x-6r_)(x + 6r_)(l + 62ri+x2)' 

The last expression is singular in the limit r ^ r_ , so F2 diverges across the inner horizon. 
The same holds for many other second-order functions, like K4 or K^. 



-26- 



References 



[1] J. Casalderrey-Solana, H. Liu, D. Matcos, K. Rajagopal, and U. A. Wiedemann, 



"Gauge/String Duality, Hot QCD and Heavy Ion Collisions," |arXiv : 1101 . 0618 [hep-th] 



[2] S. Bliattacharyya, V. E. Hubeny, S. Minwalla, and M. Rangamani, "Nonlinear Fluid 
Dynamics from Gravity," \JHEP 02 (2008) 045| , ^rXiv : 0712 . 2456 [hep-th] 



[3] R. A. Janik and R. B. Peschanski, "Asymptotic perfect fluid dynamics as a consequence of 
AdS/GFT," \Phys. Rev. D73 (2006) 045013| , |arXiv:hep-th/051216j . 

[4] N. Banerjee, J. Bhattacharya, S. Bhattacharyya, S. Dutta, R. Loganayagam, et ai, 



■'Hydrodynamics from charged black branes," JHEP 1101 (2011) 094, 



arXiv: 0809. 2596 



[hep-th] 



[5] J. Erdmenger, M. Haack, M. Kaminski, and A. Yarom, "Fluid dynamics of R-charged black 



holes," JHEPOmi (2009) 055, arXiv: 0809 .2488 [hep-th] 



[6] S. Bhattacharyya et ai, "Forced Fluid Dynamics from Gravity," JHEP 02 (2009) 018 



arXiv: 0806. 0006 [hep-th] 



[7] S. Bhattacharyya, R. Loganayagam, L Mandal, S. Minwalla, and A. Sharma, "Conformal 



Nonlinear Fluid Dynamics from Gravity in Arbitrary Dimensions," JHEP 12 (2008) HE, 



arXiv: 0809 .4272 [hep-th] 



[8] A. Chamblin, R. Emparan, C. V. Johnson, and R. C. Myers, "Charged AdS black holes and 
catastrophic holography," \Phys.Rev. D60 (1999) 064018| , |arXiv:hep-th/9902170 [hep-th] 



[9] D. T. Son and P. Surowka, "Hydrodynamics with Triangle Anomalies," Phys. Rev. Lett. 103 



(2009) 1916011 , larXiv: 0906 .5044 [hep-th] 



[10] R. Loganayagam, "Entropy Current in Conformal Hydrodynamics," JHEP (2008) 087 , 



arXiv : 080 1 . 370 1 [hep-th] 



[11] S. Bhattacharyya et ai, "Local Fluid Dynamical Entropy from Gravity," JHEP 06 (2008) 



055, larXiv: 0803. 2526 [hep-th] 



[12] P. Romatschke, "Relativistic Viscous Fluid Dynamics and Non-Equilibrium Entropy," Class 



Quant. Grav. 27 (2010) 025006, arXiv : 0906 . 4787 [hep-th] 



[13] K. Jensen, M. Kaminski, P. Kovtun, R. Meyer, A. Ritz, et at, "Towards hydrodynamics 



without an entropy current," \Phys.Rev.Lett. 109 (2012) 101601 



arXiv: 1203.3556 



[hep-th] 



[14] N. Banerjee, J. Bhattacharya, S. Bhattacharyya, S. Jain, S. Minwalla, et at, "Constraints on 
Fluid Dynamics from Equilibrium Partition Functions," \JHEP 1209 (2012) 04^ 



arXiv: 1203.3544 [hep-th] 



-27- 



[15] I. Booth, M. P. Heller, and M. Spalinski, "Black brane entropy and hydrodynamics: the 



boost-invariant case," Phys. Rev. D80 (2009) 126013, arXiv : 0910 . 0748 [hep-th] 



[16] I. Booth, M. P. Heller, and M. Spalinski, "Black Brane Entropy and Hydrodynamics," Phys 



Rev. D83 (2011) 061901, arXiv : 1010 . 6301 [hep-th] 



[17] I. Booth, M. P. Heller, G. Plewa, and M. Spalinski, "On the apparent horizon in fluid-gravity 



duality," Phys. Rev. D83 (2011) 106005, arXiv: 1102.2885 [hep-th] 



[18] G. Plewa and M. Spalinski, "In preparation," . 

[19] R. C. Myers, M. F. Paulos, and A. Sinha, "Holographic Hydrodynamics with a Chemical 



Potential," JHEP 06 (2009) 006, arXiv: 0903. 2834 [hep-th] 



[20] E. Witten, "Anti-de Sitter space, thermal phase transition, and confinement in gauge 
theories," Adv. Theor. Math. Phys. 2 (1998) 505-532, |arXiv : hep-th/980313l| . 

[21] M. Henningson and K. Skenderis, "The holographic Weyl anomaly," JHEP 07 (1998) 023, 



arXiv : hep-th/9806087 



[22] V. Balasubramanian and P. Kraus, "A stress tensor for anti-de Sitter gravity," Comrnun. 



Math. Phys. 208 (1999) 413-428| , |arXiv:hep-th/9902121 



[23] S. S. Gubser and I. Mitra, "Instability of charged black holes in Anti-de Sitter space," 



arXiv:hep-th/0009126 [hep-th] 



[24] S. Bhattacharyya, S. Lahiri, R. Loganayagam, and S. Minwalla, "Large rotating AdS black 



holes from fluid mechanics," \JHEP 09 (2008) 054 |arXiv: 0708. 1770 [hep-th] 



[25] D. E. Kharzeev and H.-U. Yee, "Anomalies and time reversal invariance in relativistic 



hydrodynamics: the second order and higher dimensional formulations," Phys. Rev. D84 
(2011) 04502^ , larXiv: 1105. 6360 [hep-th] . 



[26] R. Loganayagam, "Anomaly Induced Transport in Arbitrary Dimensions," arXiv : 1106 . 0277 



[hep-th] 



[27] K. Peeters, "A field-theory motivated approach to symbolic computer algebra," CoRR 
abs/cs/0608005 (2006) . 

[28] K. Peeters, "Introducing Cadabra: A symbolic computer algebra system for field theory 



problems," arXiv:hep-th/0701238. 



-28- 



